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Introduction.  Recent  considerations  of  the  problem  of  the  end-loaded 
cantilever  beam,  with  the  conditions  of  loading  prescribed  in  terms  of  dis- 
placements rather  than  in  terms  of  stresses,  have  lead  to  defining  relations 
for  shear  center  and  twist  center  location  in  terms  of  influence  coefficients, 
of  a particularly  simple  nature  [4,  5].  It  was  furthermore  shown  that  approxi- 
mate values  of  these  influence  coefficients,  leading  to  approximate  expres- 
sions for  shear  and  twist  center  coordinates,  could  be  obtained  by  using  St. 
Venant-type  torsional  and  flexural  stress  distributions  in  a Rayleigh-Ritz 
sense  in  conjunction  with  the  principle  of  minimum  complementary  energy  [4|. 
In  what  follows  we  extend  these  results  in  several  directions. 

We  begin  by  making  explicit  the  distinction  between  expressions  for 
shear  and  twist  center  coordinates  in  terms  of  flexibility  coefficients  (which 
were  previously  considered)  and  in  terms  of  stiffness  coefficients  (where  it 
is  shown  that  for  the  case  of  principal  centroidal  coordinate  axes  the  final 
formulas  are  as  simple  as  the  formulas  in  terms  of  flexibility  coefficients). 

We  supplement  our  earlier  statement  of  a minimum  complementary 
energy  equation  for  the  case  of  prescribed  rigid-body  type  in-plane  end 
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section  displacements  (4)  by  a statement  of  the  associated  minimum  potential 
energy  equation.  Furthermore,  we  consider,  in  addition  to  the  case  of  prescribed 
in-plane  end  displacements,  conditions  of  loading  in  a somewhat  unconven- 
tional fashion,  specifying  the  form  but  not  the  magnitude  of  these  displacements 
at  tlie  loaded  end  of  the  beam,  and  at  the  same  time  specifying  the  resultants 
but  not  the  distribution  of  in-plane  stresses  at  this  end.  We  find  that  it  is  a 
simple  matter  to  state  a minimum  potential  energy  equation  in  such  a way  as 
to  apply  to  the  latter  case,  but  we  leave  open  the  question  of  an  appropriate 
minimum  complementary  energy  equation. 

We  use  minimum  potential  and  complementary  energy  equation  statements 
for  the  establishment  of  upper  and  lower  bound  relations  for  quadratic  forms 
involving  flexibility  and  stiffness  coefficients,  respectively.  We  give  both 
types  of  bounds  for  stiffness  coefficients,  but  a lower  bound  only  for  flexibility 
coefficients,  pending  formulation  of  a minimum  complementary  energy  equation 
for  the  mixed  displacement-stress  boundary  condition  case  described  above. 

In  applying  the  principle  of  minimum  potential  energy  for  the  approxi- 
mate determination  of  flexibility  and  stiffness  coefficients  we  utilize  dis- 
placement approximations  which  have  previously  been  used  for  the  analysis 
of  the  problems  of  torsion  and  flexure  with  end-section  restraint  against 
warping  [l,7,8|.  We  think  that  it  has  not  previously  been  recognized  that 
these  approximations  may  be  utilized,  beyond  allowing  assessments  of  the 
effect  of  end  section  restraint,  for  the  purpose  of  deducing  bound  relations 
for  coefficients  entering  into  the  formulas  for  shear  and  twist  center 


coordinates. 


In  applying  the  principle  of  minimum  complementary  energy  we  signifi- 
cantly modify  our  earlier  procedure  [4)  by  not  departing  from  a St.  Venant- 
type  shear  distribution.  Instead  we  utilize  an  approximation  in  which  the  dis- 
tributions of  both  shear  and  normal  stress  over  the  cross  section  are  deter- 
mined through  use  of  the  variational  equation.  In  so  doing  the  problem  is, 
without  special  effort,  solved  for  a more  general  class  of  cases  than  hereto- 
fore, in  a manner  which  is  thought  to  represent  a significant  simplification 
of  the  earlier  work. 

Remarkably,  the  approximate  results  for  flexibility  and  stiffness  coef- 
ficients which  are  obtained  on  the  basis  of  making  quite  disimilar  approximate 
assumptions  in  connection  with  the  use  of  the  minimum  potential  energy  equa- 
tion and  of  the  minimum  complementary  energy  equation  lead  to  identical 
approximate  expressions  for  the  coordinates  of  the  center  of  shear  and  of 
twist.  This,  in  conjunction  with  our  upper  and  lower  bound  relations,  leads 
to  the  conclusion  that  our  approximate  results  are  in  fact  exact,  in  the  limit 
of  vanishing  a/L.  where  a is  a representative  cross  sectional  width  dimension 
and  L is  the  axial  length  of  the  beam,  with  the  possibility  left  open  to  refine 
the  analysis  so  as  to  account  for  the  (generally  small)  effect  of  finite  values 
of  a/L. 

A Formulation  of  the  Problems  of  Torsion  and  Flexure.  We  consider  a 
body  with  boundaries  defined  by  a cylindrical  surface  f(x,y)  = 0 and  two  planes 
z = 0 and  z = L.  We  designate  displacements  by  u,  v,  w and  stresses  by 
0x,  T , etc.  and  we  assume  that  the  normal  three-dimensional  homogeneous 


equations  of  linear  elasticity  hold.  We  further  assume  that  the  boundary 


portion  f ■=  0 is  traction  free  and  that  the  boundary  portion  z = 0 is  fixed. 

In  regard  to  the  boundary  portion  z = L we  assume  the  absence  of 

/ 

normal  tractions,  in  conjunction  with  a rigid  body  translation  and  rotation 
distribution  of  tangential  displacements,  that  is,  we  stipulate  the  conditions 

/,  = L;  a = 0,  u - IJ  - y©,  v = V + x©  . (1) 

z 

We  now  observe  that  in  writing  equation  (1)  we  may,  or  we  may  not,  stipulate 
additionally  the  magnitudes  of  U,  V,  ©.  If  we  do,  as  we  have  done  earlier 
in  conjunction  with  applications  of  the  principle  of  minimum  complementary 
energy  [4,5],  then  equations  (1)  are  a complete  statement  of  loading  conditions. 

If  we  do  not  and  leave  the  magnitudes  of  U,  V,  ©unspecified  then  we  must, 
in  order  to  complete  the  statement  of  loading  conditions  prescribe  additionally 
the  magnitude  of  two  transverse  force  components  P,  Q,  and  of  an  axial 
torque  T,  as  follows 

z=L;  f(Tx*’  vdS  = (P,Q)'  ,r(vx*  v)ds  = T • (2) 

In  view  of  the  linearity  and  homogeneity  of  the  problem  we  have  that  P, 

Q and  T will  be  combinations  of  U,  V and  6,  and  vice  versa,  in  the  form 


p = kpuu  f kpvv  + KP©e- 


Q * KQUU  + 


t = ktuu  + ...  +KTee, 


and 


upp  + cuqq  + cutt’ 


PP  + 


(3) 


(4) 
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with  the  Kpy,  etc.  being  stiffness^  coefficients,  the  Cyp,  etc.  being  flexibilit 
coefficients,  and  with  the  expectation  of  symmetry  for  the  matrix  of  the  K's 
as  well  as  for  the  matrix  of  the  C's. 

Having  equations  (3)  and  (4)  we  obtain  the  coordinates  of  the  center  of 
twist,  xT,  y , as  the  coordinates  *,y  of  that  point  in  the  end  cross  section 
for  which  u = v = 0 in  equation  (1)  while  at  the  same  time  P = Q = 0,  that  is, 
in  the  form 

yT  = XT  = ’ *V/®*P=Q  = 0 • 

and  we  obtain  the  coordinates  of  the  center  of  shear,  x„,  y_,  as  the  coordi- 
nates  of  the  point  of  intersection  of  the  lines  of  action  of  the  forces  P,  Q for 
the  case  ©=  0,  in  conjunction  with  the  torque  T being  solely  due  to  the  forces 
P,  Q,  that  is,  upon  setting  in  equations  (3)  or  (4) 

©=  0,  T = Oxs  - P yg  • (6) 

It  turns  out  that  with  these  defining  relations  the  simpler  form  of  the 
results  appears  through  use  of  the  flexibility  coefficients,  namely 

y y x ..fXT  . 

yT  C_  ’ YS  C_ _ ’ T C __  ’ S C_._  ’ 

©T  ©T  ©T  ©T 

with  y = y„  and  x = x for  the  normal  case  of  a symmetric  C-matrix. 

1 O I O 

The  corresponding  relations  in  terms  of  the  coefficients  K in  equation  (3) 
come  out  to  be  ratios  of  certain  second  order  minors  of  the  third  order 
determinant  of  the  coefficient  matrix  in  (3).  It  will  be  useful  to  note  for 


what  follows  that  we  have,  on  the  basis  of  equations  (5)  and  (3), 


KpvKoe'  kqvkp©  kpukq©'  kqukp0 

Vt  ""  kpukqv  - kqukpv  * 


1 kpukqv  ' KovKpv 


again  with  y = y and  x = x in  normal  circumstances,  and  with  the  im- 
S * Si 

portant  special -rase  formulas 


Kpe  _ Koe 

xt  r K, 


K 


(9) 


P1I 


QV 


which  result  upon  assuming  that  = ^QU  = 

Minimum  Complementary  and  Potential  Energy  Equations  for  the 

. . _ — — — — — — ■ — — _ — — — — — — — — . 


Problem  of  Torsion  and  B'lexure.  We  now  assume  that  the  material  of  the 


beam  is  such  that  its  stress  strain  relations  may  be  written  in  the  alterna- 
tive forms 


o = dA/de  , t = dA/dy  

x x xy  xy 


(10a) 


and 


€ = dB /ba  , y = SB/^t  ..... 
x x xy  xy 


(10b) 


We  then  have  for  the  case  of  prescribed  U=U,  V=V,  ©=i?as  minimum 
complementary  energy  condition  the  variational  equation  61^  - 0,  where 


I = - J J BdSdz  + UPtVQt?T. 


(11a) 


In  (Ha)  the  stresses  o , t , etc.  must  satisfy  the  differential  equations  of 
x xy 


equilibrium  and  all  stress  boundary  conditions,  and  the  variational  equation 

-f>- 


is  equivalent  to  all  strain  displacement  relations  and  displacement  boundary 
conditions.  We  have  earlier  considered  the  application  of  this  variational 
problem  for  the  approximate  determination  of  flexibility  coefficients  in  con- 
junction with  stress  distributions  corresponding  to  the  solutions  of  the  St. 
Venant  torsion  and  flexure  problem  [4,  5], 

For  an  alternate  formulation  of  Hie  problem,  within  the  context  of  the 
principle  of  minimum  potential  energy,  which  we  have  not  considered  pre- 
viously, we  now  prescribe  P = P,  Q - Q,  T = T,  in  association  with  the  end 
displacement  distribution  (1).  We  then  have  that  the  appropriate  form  of  the 

principle  of  minimum  potential  energy  is  the  variational  equation  61  , = 0 

d 

where 

Id  = JJ AdSdz  - PU  - QV  - T0.  (lib) 

In  (lib)  the  strains  e^,  y , etc.  are  given  in  terms  of  displacement  deriva- 
tives, the  displacement  components  u,  v,  w must  vanish  for  z = 0,  and  u,  v 
must  be  as  in  equation  (1)  for  z = L,  with  no  restrictions  imposed  on  U,  V.  0. 
The  variational  problem  as  stated  has  as  Euler  equations  the  differential 
equations  of  equilibrium  in  the  interior  and  all  conditions  of  prescribed  stress 
on  the  surface. 

Having  previously  used  the  variational  equation  61^  = 0 for  the  deter- 
mination of  flexibility  coefficient  approximations  - without  regard  to  the  fact 
that  the  form  of  equations  (11a)  and  (3)  indicates  that  it  would  be  more  natural 
to  use  this  relation  for  the  determination  of  stiffness  coefficients  - we  have 
not  previously  used  the  relation  61^  = 0 which,  it  is  apparent  from  (lib)  and 
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14),  is  a natural  starting  point  for  the  approximate  determination  of  flexibility 
c ocffic  icnts . 

Previous  general  considerations  on  upper  and  lower  bound- determina- 
tions for  inf'  jnce  coefficients  [2]  indicate  that  the  use  of  I is  associated  with 

s 

the  possibility  of  determining  lower  bounds  for  stiffness  coefficients  K and 
that  the  use  of  I gives  the  possibility  of  determining  lower  bounds  for  flexibility 
coefficients  C. 

4c 

Furthermore,  we  know  that  the  use  of  a potential  energy  function  I 

Cl 

defined  by 

I(,  J"J  AdSd*,  (12) 

with  I * differing  from  1^  by  the  stipulation  that  in  it  U = 0,  V = V,  0 = 0,  will 
be  involved  in  the  determination  of  upper  bounds  for  stiffness  coefficients. 

In  order  to  obtain  upper  bounds  for  flexibility  coefficients  we  should 
have  a counterpart  I to  I as  defined  in  (lla),  with  U,  V,  0 replaced  by 

B 3 

P,  Q,  T,  in  such  a way  that  the  form  of  the  tangential  end  displacement  dis- 
tribution remains  prescribed  in  accordance  with  equation  (1).  We  do  not, 
at  this  time,  know  the  way  in  which  to  introduce  these  "partial"  displacement 

boundary  conditions  into  the  principle  of  minimum  complementary  energy. 

♦ 

Because  of  this  we  do  not  here  establish  the  appropriate  form  of  I (which, 
if  symmetry  considerations  were  the  principal  guides,  ought  to  be  given  by 

- JjBdSdz). 

Upper  and  Lower  Bound  Relations  for  Influence  Coefficients.  Appro- 
priate transformations  of  the  functionals  1 , I . and  I . , as  defined  In 

8 a d 
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equations  (11)  and  (12),  lead  to  the  upper  and  lower  bound  relations 


I * £(UP  + VQ  t ?T)  s r * 
8 d 


(13) 


- i(up  + vq  + er>  s i , 

Cl 


(14) 


with  the  missing  left  hand  side  of  equation  (14)  making  it  evident  that  it  could 

* 

be  useful  to  have  a functional  1^  as  discussed  at  the  end  of  the  preceding 
section.  In  equation  (13)  I corresponds  to  the  functional  I in  (lla),  with 

8 S 

any  stresses  a^,  T , ctc>  which  satisfy  equilibrium  differential  equations  and 

~ * * 

stress  boundary  conditions,  and  1 corresponds  to  I in  (lib),  with  any  d if  - 

d d 

ferentiable  displacement  state  u,  v,  w which  satisfies  the  stipulated  displace- 
ment boundary  conditions,  with  the  same  rules  connecting  I in  (14)  and  I 

d d 

in  (lib). 

In  order  to  see  that  (13)  and  (14)  represent  bound  relations  for  stiffness 
and  flexibility  coefficients  respectively,  we  observe  that  I ^ as  well  as  1 ^ will 
be  quadratic  forms  in  U,  V and  ©,  which  may  be  written  as 


l*  ^KtUtItU2  + K”vUVf... 


uu 


ee 


(15a) 


^Ke©e2 


(15b) 


At  the  same  time  the  quantity  which  is  bounded  from  above  and  below  in  equa- 
tion (13)  may  be  written,  with  the  help  of  equations  (3),  ns 

UP  * VQ  ♦ eT  . Kpuu2  t (Kpv  < kqu>uv  I ...  t Krpe2  llfcl' 
KPU  ■ K„u-  *IKPV  * kqu’  ■ KPV  ■ KQU  * Kuv  »•  «“« 


The  author  would  like  to  acknowledge  that  the  result  expressed  by  equations 
(13),  (15)  and  (16)  was  found  independently  by  S.  Nair. 
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when  V - & = 0 equations  ( I J) , (I1’)  and  ( 1 (>)  specialize  to  a previous]  y slater) 
result  for  the  problem  of  flexure,  within  the  framework  of  the  theory  of  plane 
stress  [3].  Analogously,  setting  U = 0 and  V = 0 gives  a simple  explicit  result 
for  the  problem  of  twisting  with  end  restraint  against  warping,  which  may  not 
have  been  stated  previously.  ^ We  further  note  the  evident  simplifications 
which  occur  in  tile  above  upon  stipulating  symmetry  conditions  which  result 
in  bound  relations  for  quadratic  forms  in  two  variables  in  place  of  the  relations 
for  the  general  three-variable  case. 

We  next  consider  equation  (14),  written  in  the  form  - I £ &(UP  + VQ  + PT). 

d 

In  this,  equation  (lib)  enables  us  to  write  with  suitable  coefficients  C^J, 


* cpqPq  f 


+ JCttt2’ 


(17) 


At  the  same  time,  we  have,  through  the  use  of  (4),  that 


pu  t qv  t re 


C pz 
UP 


* <°UQ  + CVP,PQ 


+ V2’ 


(18) 


with  the  obvious  consequences  that 


°PP  * CUP’  C'QQ  * CVQ’  CTT  * CeT* 


(19) 


but  with  the  determination  of  bounds  for  the  coefficients  + Cyp,  etc. 

of  the  mixed  terms,  and  therewith  of  bounds  for  the  coordinates  of  the  centers 
of  shear  and  of  twist,  depending  upon  the  establishment  of  a bound  functional 


I for  the  lefthand  side  of  equation  (14). 


The  result  for  twisting  without  restraint  against  warping  is  contained  in  an 
early  fundamental  paper  by  Trefftz  [6|. 
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Approximate  Determination  of  Shear  and  Twist  Center  Location  Through 
Use  of  Minimum  Potential  Energy  Principle.  Similar  to  what  i as  been  done 
in  earlier  work  on  problems  of  combined  twisting  and  bending  of  beams  [1,  7, 
8]  we  begin  by  stipulating  as  approximations  for  components  of  displacement 

u = u(z)  - y8(z),  v - v(z)  f x9(zl  , (20) 


w = w (z)  ♦ xar(z)  ( yfl(z)  + 0(x,  y)X(z)  (21) 

o 

In  this  0(x,y)  is  a function  which  is  to  be  assumed  suitably,  with  the  various 
functions  of  z in  (20)  and  (21)  to  be  determined  by  the  variational  procedure, 
with  or  without  imposition  of  additional  constraint  relations. 

As  regards  stress  strain  relations  we  here  consider  a material  pos- 
sessing a limiting  type  orthotropy,  in  such  a way  that  there  is  partial 
rigidity,  with  each  cross  section  translating  and  rotating  as  an  entity.  * 

Our  limiting-type  stress  strain  relations  are  three  relations  i - ( -y 

x y xy 

= 0,  in  conjunction  with  three  relations  of  the  form 


a 

z 


= Ec  , 

z 


T 

yz 


Gy 


Y7‘ 


(22) 


With  (22)  we  have  as  expression  for  the  strain  energy  function  A in  equation 

(ID 


(23) 


This  assumption  is  meaningful  for  sufficiently  slender  beams  only,  where 
its  approximate  validity  depends  on  the  relative  insignificance  of  the  com- 
ponents of  stress  0*.  tTy,  TXy  in  comparison  with  the  components  oz,  txz. 
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where,  on  the  basis  of  (20)  and  (211, 


e = w ' v xa'  4 yP'  4 0X',  (24) 

z o 

y = u'  + a - y9’  4 X#  , Y = v'  4 $ + x0'  + X0  . (25) 

xz  ,x  yz  ,y 

In  introducing  (23)  to  (25)  into  the  strain  energy  integral  we  shall  assume 
that  the  origin  of  the  x,  y system  of  axes  is  at  the  elastic  centroid  of  the  cross 
section,  that  is  we  stipulate  the  relations  j"(x,  y)EdS  = 0.  Furthermore,  we 
shall  assume  that  0 is  the  warping  function  for  St.  Venant  torsion  of  a homo- 
geneous beam  with  the  same  cross  section  as  the  given  beam,  that  is,  we 

i shall  assume  that  0 is  determined  through  the  relations  V^0  = 0,  f(0  - y,  0 

, x , y 

+ x)dS  = (0,  0),  (-y  t 0 )dy  - (x  + 0 )dx  * 0 along  the  boundary  f = 0 of  the 

» x » y 

cross  section,  we  shall  set  as  an  abbreviation 


V = ^(0  2 + 0 2)dS  = J (y 0 - x0  )dS  (26a) 

» x f y *y 

and  we  shall  stipulate,  as  we  may,  that  J'E0dS  = 0.  To  be  consistent  with  our 

choice  of  0 we  furthermore  assume  in  what  follows  that  G - const.  With  this 

we  now  obtain  as  expression  for  the  approximation  1^  to  the  functional 

I in  equation  (11) 
d 

T.  = §1  {S  (w  ')l  4 I (a')2  4 I <P')2  + 21  o'P' 
d 1 E o xx  yy  xy 

+ nx')2  + 2r  a'\'  4 2r  p'x' 

x y 

4 G[(u'  4 o)2  4 (v'  4 P)Z  - 2(u'  4 or)9' /ydS 
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+ 2(v'  + P)0',fxdS  + I (0')2  + DX2 
p 

- 2DX0']}dz  - PU  - QV  - T©  , (27) 

where  U = u(L),  V = v(L),  © = 0(L),  and 

i =J(xZ  + y2)ds  . <r,r  T)  = |(i,x,y)tos  , 

p x y 

(SE’Ixx’Ixy,Iyy)  = XV’  y2)EdS  * <26b) 

Inasmuch  as  we  are  concerned  with  approximate  rather  than  exact 
results,  we  shall  now  further  assume  that  translational  deflections  due  to  trans- 
verse shear  may  be  neglected  and  that  the  entire  transverse  shear  strain  energy 
is  that  due  to  twisting.  Considering  the  form  of  (27)  the  desired  reduction  will 
be  accomplished  upon  introducing  the  additional  constraint  relations 

or  = — u , 0 = -v'  , X = 6'  . (2fi) 


Anticipating  furthermore  the  result  w^  - 0,  we  will  then  have  in  place  of  equation 


(27) 


I . =7J*fl  (u')2  + I (v*)2  + 21  u'v" 
d 2 • xx  ' yy  xy 


+ r<ev  - zr  e'u*  - 2r  0%* 

x y 


+ C(0'r)dz  - PU  - QV  - T©  , 


(29) 


with  C = 0(1^  - D)  being  the  conventional  St.  Venant  torsional  stiffness  factor. 

In  evaluating  the  variational  equation  6l  = 0 we  take  account  of  the 

d 

constraint  boundary  conditions  u(0)  = v(0)  = 0(0)  = u'(0)  = v'(0)  = B' {0)  = 0. 

The  corresponding  six  conditions  for  the  loaded  end  of  the  beams  are  the  Euler 
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boundary  conditions 


i i i - r e*  - o . i + 1 v"  - r 0"  = -p 

xx  xy  x xx  xy  x 


I u'ti  v*  - r 0*  = 0 , I u"  4 I v" - r 0"  = -Q 
xy  yy  y xy  yy  y 


(30) 


r .1'  4 r v'  1 re*  = 0 , -r  . r v"  4 re"  -c0'  = -t 

x y x y 

for  7.  - L,  with  the  Euler  differential  equations  being 

t IV  IV  r ftiv 

Iu  4 I v - T W =0  , 

xx  xy  x 


, lV  , r IV  r flIV  n 
I u 4lv  -la  =0 


xy 


yy 


(31) 


-r  UIV  - r vIV  4 r^v  - c0*  = 0 . 

x y 

In  order  to  solve  the  problem  as  stated,  we  begin  by  obtaining  from 
(31)  the  transformed  differential  equations 


Kurv  = <r  1 . r 1 »0IV,  kviv  = <r  1 - r 1 >eIV 

x yy  y xy  y xx  x xy 


r*0iv-c0"=o  • 


(32) 

(33) 


where 


k = 1 i -i2  . r =r-<r2i  -zrn  4rzi  )/k  , 

xx  yy  xy  * x yy  y x xy  x xx 

and  from  (30)  the  transformed  constraint  boundary  conditions 


(34) 


u'  = o , ku"=  -i  p 4 1 q 4 (i  r - 1 r je* 

yy  xy  yy  x xy  y 


v*  = o , Kv"  = -i  q 4 1 p 4 (i  r . i r )0* 

xx  xy  xx  y xy  x 


(35) 
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a*  ^ o , r 8"  - - -t  - (i  r -i  r )(p/K)  - (i  r -i  r hq/ki  , (36> 

* yy  x xy  y ' xx  y xy  x ' ' ' 

for  7 - L. 

Equations  (13)  ;>nd  (10),  in  conjunction  with  the  conditions  0(0)  = 0 (0)  = 0, 
give  as  expression  for  0, 


0(z) 


r i - r i 

* yy njii  p , 

2 

i r - i 
xx  yy  xy 


r i - 

y xx 

I I 
xx  xy 


X 


7 

L 


sinhXL  - sinhX(L  - 
XL,  cos  hXL, 


(37) 


Having  0 as  in  (17),  we  find  u and  v from  equations  (32)  and  (35)  in  the 

form 


u(z) 


v(z) 


PI  - QI 
—11 

II  - I 
xx  yy 


QI  - PI 

xx 

II  - I 

XX  yy 


HI 

2 

xy 

HI 

2 

xy 


r i . r i 

_2L1X 5^0(7) 

II  - I 
xx  yy  xy 

r i - r i 

.XXX — i_SLe(z) 

ii  - 1 
xx  yy  xy 


(38) 


We  use  equations  (17)  and  (18)  in  particular  in  order  to  obtain  the  values 
of  & - 0(1,),  U = u(L),  V - v(L)  in  terms  of  P,  Q,  T,  so  as  to  obtain  from  equa- 
tions (24),  as  approximate  expressions  for  flexibility  coefficients 


L f,  tanhXl,\  „ L *x*yy  *Vxy  /.  tanhXL\ 

C«T  c V ' XI,  ) cep  c t x r2  ( XL  ) 

xx  yy  xy  ' ' 


» 


k 


b 


'UP 


L_ 

3 


JOL 


II  - I 

xx  yy  xy 


[xIyy.-..  ry’.xy  1 

\ i i - i2 

' xx  yy  xy 


1 - 


ta  n h XL 

XI,  J 


(39) 


etc.  Introduction  of  these  expressions  into  the  defining  relations  (7)  for  the 
coordinates  of  the  center  of  shear  and  of  twist  then  give  as  approximations 
for  these  coordinates 


ys  =yT 


r i - r i 
x yy usx 

ii  - iz 
xx  yy  xy 


r i - r i 

y xx  x xy 

II  - I2 
xx  yy  xy 


(40) 


We  note  that  in  the  event  that  the  x,  y-axes  in  the  cross  section  are 

principal  elastic  axes  we  have  I =0,  and  equations  (40)  reduce  to  the  simplified 

xy 

form  x = -T  /I  , y -T  / 1 . For  the  case  of  a constant  modulus  E,  these 

s y yy  » * xx 

latter  formulas  agree  with  the  results  previously  obtained  in  the  classical 
Weber-Trefftz  considerations,  as  well  as  with  our  earlier  approximate  results 
which  followed  from  a determination  of  flexibility  coefficients  through  use  of 
the  principle  of  minimum  complementary  energy  in  conjunction  with  approxi- 
mations for  stresses  as  given  by  the  St.  Venant  theory  of  torsion  and  flexure  1 4 j. 

An  interesting  special  case  of  the  above  is  the  case  of  a solid  circular 
cross  section  for  which  4)  ■=  0,  throughout.  We  then  have  that  the  location  of  the 
shear  center  coincides  with  the  elastic  centroid  of  the  cross  section,  for  all 
possible  variations  of  E,  as  long  as  it  is  assumed  that  G does  not  vary.  This 
result  should  be  compared  with  the  well-known  result  for  a uniform  semi- 
circular cross  section  which  may  be  interpreted  as  the  case  of  a complete 
circular  cross  section  with  vanishing  F.  and  G over  one-half  of  the  section. 
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For  this  case  we  have  as  the  distance  of  the  centroid  from  the  straight  portion 


» 


of  the  cross  sectional  boundary  curve  x = 4a/3ff  **  0.42a,  while  at  the  same 

c 

time  the  distance  of  the  shear  center  is  given  by  = 8a/5ff  **  0. 51a,  with  the 
differences  in  assumptions  concerning  the  distribution  of  G for  the  two  cases 
evidently  being  responsible  for  a significant  effect  on  the  location  of  the  center 
of  shear  and  of  twist. 


Approximate  Determination  of  Stiffness  Coefficients  Through  Use  of 

Minimum  Complementary  Energy  Principle.  We  now  consider  the  use  of  the 

variational  equation  - 0 with  1^  given  in  equation  (11a),  with  = 

T =0,  and  with  complementary  energy  density 
xy 


D - a(oI 2/ E + T2/G  + T2/G) 
x y 


(41) 


where  0 - a , T = T , T - T 

7.x  xz  y yz 

We  have  previously  considered  the  application  of  this  principle,  for 
the  case  of  cross  sections  with  the  x-axis  an  axis  of  symmetry,  and  with  E 
independent  of  x and  y (4),  on  the  basis  of  stipulating  a St.  Venant  distribution 


I tJ  r F.Q(L  - z)y  , 

yy 


(42) 


\ 'QIW  * T,.y  • Ty  *Q|,y  - *.x>  - T\,  <«> 

with  X and  ♦ as  stress  functions,  with  t and  t as  particular  solutions  of 

x y 

I (t  + t ) - Ey  = 0 , (44) 

yy  x,  x y,  y 

and  with  X , ♦ , Q and  T to  be  determined  by  the  variational  equation  in  conjunc- 


tion with  the  boundary  condition  T^dy  - T^dx  = 0 along  f(x,  y)  3 0. 

It  was  found  in  |4|,  with  the  help  of  transformations  of  some  complexity. 
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that  the  result  so  obtained,  comes  out  to  be  x = x = -T  /I  an<3  an  analogous 

91  y yy 

outcome  may  be  anticipated  for  the  problem  without  an  axis  of  symmetry,  as 
long  as  the  analysis  is  restricted  by  the  assumption  that  E is  independent  of 
x and  y. 

In  what  follows  we  consider  the  derivation  of  a different  approximate 
result,  of  independent  interest  and  — insofar  as  our  subsequent  bound  calcula- 
tions are  concerned  — associated  with  slightly  better  results  than  would  follow 
from  the  use  of  an  equation  equivalent  to  (42). 

Our  starting  assumption  is  now,  in  place  of  equation  (42),  the  relation 

where,  as  in  (43),  T and  T are  independent  of  z and  subject  to  the  constraint 
x y 

boundary  condition 


f(x,  y)  = 0 ; r dy  - T dx  ; 0 . 

x y 


(46) 


An  introduction  of  (45)  and  (46)  into  (41)  and  (11a)  leaves,  upon  carrying 
out  the  integration  with  respect  to  z,  the  approximate  energy  expression 

r = XJ[(U  - yS)Tx  + (V  + x0)Ty 
. 3 


(T  +T  )2  - 7p(TZ  + TZ)ldxdy  . 
6E ' x,x  y,  y 2G  x y 


(47) 


The  variational  equation  6lg  = 0,  with  the  stipulation  that  fir^dy  - ®T^dx  ' 0 for 
f = 0,  has  as  Euler  equations  two  differential  equations  for  and  T , of  the  form 


T 

x 

G 


LZ  fTx,x  * Ty,y\  _ U 
3\  E ) x-  L 


-IK- 


r 2 it  +r  \ — 

JL  L.  ( x,x  y.  y \ _ v_L 

G 3 \ E ) ~ L 


X0 


(48) 


Equations  (48)  may  be  simplifier!  by  deducing  from  them  the  relation 


ft)  ■('#) 

' '.y  ' x 


-z* 

L 


(49) 


which,  in  turn,  implies  as  expressions  for  T and  T in  terms  of  an  arbitrary 

x y 

function  0(x,  y). 


T =G(0  -Oy/L)  , T =G(*  +ex/L)  . 

x * y » y 

Introduction  of  equations  (50)  into  (48)  then  gives  further 


(50) 


2 /(G*  J v ♦ <G*  J 


0 - k_(. 

,X  3 \ E 


2/(G0  ) f (G0  ) 

0 . k_  i2L^ -JUJL 

.y  3\  e 


(51) 


y '.y 

While  it  is  possible  to  continue  the  analysis  for  variable  G,  the  results 
of  not  doing  tins  become  sufficiently  simpler  to  justify  a restriction  from  here  on 
to  the  case  G = const.  With  this  restriction,  and  with  observation  of  the  condi- 
tion jE0dS  - 0,  we  readily  obtain  from  (51),  as  a second  order  differential 
equation  for  0, 


0 - (GL2/3E)720  = Ux/L  + Vy/L  . 


(52) 


Equation  (52)  differs  significantly  from  the  corresponding  equation  for 
the  theory  of  torsion  and  flexure  in  accordance  with  St.  Venant  by  the  appearance 
of  the  first  term  on  the  left.  We  note  that  for  slender  beams,  with  repreaenta- 
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five  cross  sectional  dimension  a «L  this  term  will  be  small  compared  to  the 

2 2 

second  term,  of  relative  order  Ea  /GL  , and  that  considering  the  form  of  the 
differential  equation,  we  may  take  account  of  this  term  by  an  iterative  procedure. 
The  physical  reason  for  the  occurrence  of  the  first  term  in  (52)  is  evidently  the 
stipulation  of  a condition  of  no  cross  sectional  warping  at  the  fixed  end  of  the 
beam.  While  this  condition  is  disregarded  in  the  St.  Venant  formulation,  it  is 
taken  account  of,  approximately,  in  the  present  approximate  solution  through 
use  of  the  principle  of  minimum  complementary  energy. 

In  order  to  obtain  approximate  expressions  for  flexibility  coefficients 
we  now  consider  the  solution  of  (52)  to  consist  of  three  parts, 

0 =O0fl/ L + (U0  /L  t ?0  /L)/G  (53) 

H u u 

with  the  functions  0O,  0 , 0 subject  to  the  differential  equations. 

“ 11  V 


V . JE 


r)<V  *„ 


’ *v*  ’ 2 (0,  x*  y)  ' 

L 


and  to  the  boundary  conditions 


f * 0 ; 0.  dy  - 0O  dx  : ydy  t xd> 


f = 0 ; (0  ,0  )dy  - (0  ,0  )dx  = 0 

UgX  v , X r U,y  v, 


With  0O,  0 and  0 determined  through  equations  (54)  and  (55),  we  then 
have  as  expressions  for  cross  sectional  forces  P,  Q and  torque  T,  on  the 
basis  of  equations  (2)  and  (50), 

PL  =GOj(0o  - y)dS  + uj 0 dS  + vj*0  dS  . (56a) 

0,  X 11,  X v,x 
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T 


t 


with  a corresponding  expression  for  Q,  and 

TL  = G0j(x0e  y - y x + x2  + y2)dS 

t-U.f(x0  - y0  )dS  + v/(x0  -y 0 )dS  . (56b) 

u, y u,x  v,y  v,x 

A comparison  of  equations  (56)  with  equations  (3)  gives  as  approximate 
expressions  for  stiffness  coefficients 

K„,  =L'IJ0  dS,K_,  =L*,j0  dS  , 

PU  u,x  PV  J v,x 

KpQ  = CL’  x * y)dS . (57) 

K = L*  Jj*(x0  -y 0 )dS 

TU  u,  y ’ u,x' 

= GL'1J(x00  y - y0g  ^ + x2  + y2)dS  . 

Stiffness  Coefficients  Expressed  in  Terms  of  Warpins  Function  for 
St,  Venant  Torsion.  It  is  convenient  to  designate  the  warping  function  for 
St.  Venant  torsion  within  the  present  context  by  0^\  with  this  function  being 
tiie  solution  of  the  boundary  value  problem. 

7?'0^O>  = 0 ; (0g^xdy  - 0g°ydx)b  = (xdx  + ydy>b  . (58) 

W e then  have* 

S<  - y)dS  = KV  + *)dS  = 0 , (59) 

and 

f(**e°}Y  * y*(e°l  * *Z  f yZ>ds  s ip  * o = c . <6o) 

* See,  for  example,  Cove's  Treatise,  Fourth  Edition  (1934),  pp.  311-313. 
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We  note,  specifically,  on  the  basis  of  equations  (59)  that  0g  ' is  associated 
with  vanishing  values  of  the  coefficients  and  K^g  and  that  therefore 

equation  (58)  represents  an  inadequate  approximation  to  the  contents  of  equa- 
tions (54)  and  (55)  insofar  as  the  determination  of  is  concerned.  We  resolve 
this  difficulty  by  considering  the  improved  approximation 


with  determined  from  the  relations 


('•«) 


GL 


We  now  obtain,  upon  observation  of  (62), 


(62) 


Kpe  ■ol'K> 


„ * t-i'V , - -^i01 

x ,y  ‘y  GL 


dS 


= -3L_3jEx^0)dS  , 


(63a) 


with  a corresponding  relation 

kq©  = -3L"3J>Ey^e°)dS  • (63b) 

and  with  Kjq  Riven,  on  the  basis  of  (60)  and  (57),  by 

KtQ  = GL*  *C  . (63c) 

In  evaluating  the  remaining  stiffness  coefficients,  we  may  use  the 


i 
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t 


approximations  t>  - 0^°*  and  <t>  - 0^,  with  the  boundary  conditions  (55b)  and 

r u u v v 

the  differential  equations  ?2(0^,0^)  = -3EL  2(x,y).  With  this  we  obtain 


K = L'1J^,0)  dS 
PU  u,x 


= L'11T[(x0(O)  ) 4 (x0(O)  ) + 3EL~2x2]dS 

u,  x , x u,  y , y 

r l''$x(0(O)  dy  - 0(o)  dx)  4 3L_3jEx2dS 

IJ  , x II , y 

■-  31/  ,J*Ex2dS  , 


(64a) 


and.  analogously, 


KPV  = Kotl  3 3>-*3jExyHS  , K = 3L'3jEy2dS  . 


3r„  2. 


(64b) 


PV  QU  ” ' " ' "QV 

Evidently,  these  results  are  such  that  the  effect  of  transverse  shear  deforma- 
tion is  not  included,  and  it  is  apparent  that  the  calculations  including  this  effect 
will  depend  on  a consideration  of  functions  0^,  0^  in  approximations  0 = 

0(O)  + 0(l),  etc.  with  ?20(,)  = -(3E/GL2)0(O).  etc. 
u u u ' u 

It  remains  now  to  evaluate  K,j,  and  K^<y»  in  such  a way  as  to  express 
these  quantitites  in  terms  of  integrals  involving  the  function  0g^>  if  possible. 
This  is  ac  complishod  as  follows.  We  now  use,  in  the  defining  relation 

K r ..-W0)  - y0(O>  )dS  , (65) 

HI  u,  y ii,  x ' 

Green's  theorem,  involving  and  0^  and  appropriate  relations  satisfied 

by  these  functions,  in  the  form 


:<0<O)  0(o)  + 0<°>  0(o)  ,ds 

B,  x u, x B,  y u,  y 


=H°’< >y  - »!°%  • >'.0,*vo>dS 


».  y 


1 
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= 3L  '2jEx0^O)dS  . (66 

In  order  to  see  that  the  left-hand  side  of  (66)  is  in  fact  what  we  wish  to  have, 
in  place  of  the  right-hand  side  in  (65),  we  now>  make  use,  in  place  of  \ of 
the  associated  torsion  stress  function  defined  by  - y and 

, y 0.x 

♦ ^ - x,  and  therewith  by  -2  and  Ofr^dy  + '1^  ^dx)  - 0. 

. x o.  y , y , x b 


We  then  have 


M°>  0<°>  + *<°>  *(°>  )dS 
0,  X u,  x 9,  y u,  y 


,J|(y  t*(lV0)  -(x  t Ms 


. y x 


x n , y 


.,jW°>  - x*(0)  )dS  + jW>  -*"V'»,dS  . (67, 

u,x  u,  y ,y  u,y  ,x  u,  y 

A second  application  of  Green's  theorem,  now  to  the  second  integral  on  the 

right,  and  observation  of  the  properties  of  shows  that  this  integral  vanishes. 

Therewith,  and  with  (66,  and  (65),  we  then  have  as  expression  for  K.r(J  in  terms 

-*r- 

Ktu  = -3L'3jE*0(eO)dS  . (68a) 

An  analogous  reduction  gives 

Ktv  = -3L~3jEy0(en)dS  . (68b) 

We  note  from  (68)  and  (63)  that  our  stiffness  coefficients  do  satisfy,  as  they 
should,  the  symmetry  relations  K^,y  = KpQ  and  = Kq©,  and  we  also  note 

that  we  have  previously  used  a similar  transformation  as  in  going  from  (65)  to 
(68a),  in  expressing  the  approximate  value  for  = x^,  in  f4)  for  cross  sections 
symmetric  about  the  x-axls  in  terms  of  St.  Venant's  torsional  warping  function. 
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Having  equations  (60),  (61),  (64)  and  (68)  for  stiffness  coefficients, 
we  now  see  from  equations  (fi)  for  the  coordinates  of  the  centers  of  twist  and 
of  shear  in  terms  of  these  coefficients,  in  conjunction  with  the  defining  rela- 
tions (26),  that  upon  identifying  0g  ^ with  the  function  0 in  (26)  the  present 
approximate  analysis  by  means  of  the  principle  of  minimum  complementary 
energy  does  in  fact  lead  to  the  same  equations  (40)  for  the  location  of  the 
centers  as  obtained  through  use  of  the  principle  of  minimum  potential  energy. 

It  is,  however,  worth  noting  in  this  connection  that  while  in  the  complementary 
energy  calculations  the  function  0g  ^ appears  as  a logical  consequence  of  the 
analysis,  it's  corresponding  appearance  in  the  potential  energy  calculations 
depends  on  a fortuitous  ad  hoc  assumption  in  the  displacement  approximation 
equation  (21),  and  that  no  agreement  between  the  two  types  of  results  would 
have  occured  if  instead  of  defining  0 in  equation  (20)  as  St.  Venant's  torsional 
warping  function  some  other  definition  had  been  used. 

Some  Explicit  Bounds  for  Influence  Coefficients.  A return  to  our 
consideration  of  bound  relations  in  equations  (17)  to  (19)  indicates  that  the 
set  of  approximate  flexibility  coefficients  C,  in  equations  (37)  to  (39),  is  in 
fact  also  a set  of  lower  hound  coefficients  C^.  In  this  connection,  we  particu- 
larly note  the  factor  I - (Al,)  1 tanhAL  in  equations  (39)  which  makes  these 
coefficients  smaller  than  they  would  be  without  this  factor. 

Furthermore,  we  may  utilize  the  analysis  in  equations  (20)  to  (38)  for 

the  purpose  of  solving  the  analogous  problem,  with  U = U,  V = V,0  = 8 as 

constraint  conditions  for  z - L,  instead  of  the  conditions  P = P,  Q = Q,  T = 

T,  and  with  I in  (lib)  replaced  by  I in  (12),  in  such  a way  that  the  values 
n d 
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of  P,  Q,  T which  occur  in  21.  = UP  + VQ  f »T  are  the  same  combinations 

d 

of  loaded-end  values  of  derivatives  of  u,  v,  9 as  occur  in  the  expressions  for 
p(  Q,  T in  equations  (30).  It  follows  from  this  that  the  upper  bound  stiffness 
coefficients  KU  in  equations  (15a)  are  in  fact  the  elements  of  a matrix  which 


is  the  inverse  of  the  matrix  C^which  is  implied  by  equations 


Having  thus  obtained  an  upper  bound  quadratic  form  for  the  coefficients 
K of  the  matrix  K,  in  accordance  with  equations  (15a)  and  (13),  we  next  observe 
that  our  analysis  in  equation  (45)  to  (57)  is  of  such  nature  as  to  make  the  ele- 
ments K in  equations  (57)  effectively  lower  bound  coefficients  K^,  in  accordance 
with  equations  (15b)  and  (13).  Evidently,  equations  (53)  to  (57)  no  more  than 
enable  us  to  calculate  these  coefficients  K^.  However,  considering  the  form 
of  the  differential  equation  (52)  and  the  order  of  magnitude  considerations  leading 
from  (54)  to  (58),  (f>l),  (62)  and  (64),  wo  have  that  the  explicit  approximations 
for  the  coefficients  K in  equations  (63),  (64)  and  (68)  are  in  fact  the  values  of 
the  lower  bound  coefficients  K^,  except  for  additive  terms  of  relative  order 
(a/L)2. 

We  may  associate  this  conclusion,  with  another  one  which  will  result 
from  a consideration  of  the  upper  bound  coefficients  K obtained  by  inversion 
of  the  matrix  C*J.  This  complementary  conclusion  is  that  the  upper  bound 
coefficients  KU  obtained  in  this  manner  also  agree  with  the  approximate  coef- 


ficients in  equations  (63),  (64)  and  (68),  except  for  terms  of  relative  order 
2 

(a/L)  . This  being  the  case  it  is  then  possible  to  state  that  the  approximate 
values  of  the  coefficients  K in  (63),  (64)  and J68)  are  in  agreement  with  the 


exact  values,  except  for  terms  of  relative  order  (a/L)  . It  follows  then  further 
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that  equations  (40)  for  Hit  coordinates  of  the  centers  of  twist  and  of  shear 
represent  the  location  of  these  centers — as  defined  by  equations  (5)  and  (6),  in 
association  with  the  described  mixed  boundary  value  problem  in  three-dimen- 
sional linear  elasticity  theory  — exactly,  except  for  terms  of  relative  order 
(a/U)^.  We  may,  if  we  wish,  obtain  improved  bounds  K^,  including  terms 
of  relative  order  (a/I.,)^,  by  extending  the  calculations  based  on  equations 
(51)  to  (57)  to  the  extent  of  determining,  by  iteration,  the  functions  0^  + 

*<"  + 0^  + 0^,  0^  + 0^.  We  may  also  obtain  improved  bounds 

by  carrying  out  the  analysis  based  on  the  displacement  approximations 
(20)  without  imposing  the  constraint  relations  (28),  In  contemplating  such  a 
program  it  must,  however,  be  born  in  mind  that  the  present  analysis  is  based 

on  the  assumption  of  a medium  for  which  the  axial  stress  a results  in  no 

z 

lateral  contraction  effects  in  x,y-planes.  Consideration  of  this  lateral  contrac- 
tion effect  would  mean  greatly  increased  complexity  of  the  calculations  leading 
to  the  values  K^.  Such  calculations  may  be  expected  to  leave  the  first-approxi- 
mation bound  results  unchanged,  while  at  the  same  time  being  responsible  for 
second-approximation  bound  results  involving  additional  terms  of  relative  order 


a/L  as  well  as  of  order  (a /L)^. 
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